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CHAPTER  I 

SUMMARY  AND  CONCLUSIONS 


This  project  report  discusses  Phase  II  of  the  work 
effort  under  Contract  DACW63-77-C-0003  with  the  Fort  Worth 
District,  U.S.  Army  Corps  of  Engineers.  This  work  effort 
involved  the  development  and  application  of  predictive 
models  to  evaluate  the  impact  of  the  proposed  Tennessee 
Colony  Lake  on  the  Trinity  River,  Texas.  Specifically, 

Phase  II  deals  with  the  interim  pool  configuration  of  the 
lake  (elevation  240  feet) ,  whereas  Phase  I  dealt  with  the 
final  configuration  (elevation  275  feet).  The  Phase  II 
effort  examined  the  predicted  water  quality  and  eutrophi¬ 
cation  conditions  in  the  lake  and  evaluated  the  predicted 
dissolved  oxygen  conditions  in  the  lake  headwaters  and 
in  the  Trinity  River  below  the  dam. 

The  following  paragraphs  summarize  the  major  results 
of  Phase  II  of  this  study: 

Model  Configurations 

1.  The  geometry  of  the  interim  pool  was  estimated  using 
available  topographic  maps  and  channel  cross-sectional 
data.  Two  configurations  were  developed  for  analysis: 
a)  a  completely-mixed,  single-segment  configuration 
with  a  total  volume  of  109,900  acre-ft,  a  total  sur¬ 
face  area  of  15,100  acres,  and  an  average  depth  of 
7.3  feet;  and  b)  a  six-segment  model  with  three  deep 
main  channel  segments  and  three  shallow  bay  segments. 
The  six  segments  consist  of  one  channel  and  two  shallow 
segments  in  the  lower  reservoir,  one  channel  and  one 
shallow  segment  in  the  middle  reservoir  and  one  upper 
channel,  or  headwater,  segment. 


2.  The  two  annual  flow  regimes  developed  in  the  Phase 
I  effort — a  low  May  flow  and  a  high  May  flow — 
were  used  as  the  bases  for  evaluating  the  interim 
pool.  Estimates  of  daily  loadings  of  BOD  ^ ,  nitrogen, 
phosphorus,  and  suspended  solids  from  the  Trinity 
River  inflow  were  based  upon  rating  curves  of  con¬ 
centrations  versus  flow  developed  in  Phase  I. 

Chlorophyll  'a'  Simulations 

3.  The  eutrophication  model  framework,  which  was  cali¬ 
brated  and  validated  in  Phase  I  for  Lake  Livingston, 
was  applied  to  the  proposed  Tennessee  Colony  Lake 
for  both  the  completely-mixed  and  six-segment  con¬ 
figurations  during  low  and  high  flow  hydrologies. 

The  models  predicted  daily  concentrations  of  nutrients 
solids,  and  chlorophyll  'a'  in  the  completely-mixed 
segment  and  each  of  the  six  separate  segments.  Chloro 
phyll  'a*  was  previously  determined  to  be  an  adequate 
indicator  of  algal  numbers  and  of  primary  productivity 
The  sensitivity  of  the  six-segment  model's  chlorophyll 
'a'  results  to  the  rate  of  horizontal  dispersion  was 
also  examined  by  making  alternate  model  runs  with 
the  dispersion  rate  reduced  to  the  lower  bound  value 
for  lakes. 

4.  The  six-segment  model  indicates  that,  under  low  May 
flow  hydrology,  chlorophyll  'a'  concentrations  in 
various  sections  of  the  reservoir  are  projected  to 
range  from  5  to  130  yg/l.  Concentrations  in  the  upper 
channel  segments  remain  at  relatively  low  levels  "ntil 
midsummer,  at  which  time  substantial  growth  to  levels 
of  approximately  100  yg/l  is  expected  to  occur.  By 
contrast,  in  the  shallow  segments  high  levels  of 
chlorophyll  'a',  between  60  and  130  yg/l,  are  pro¬ 
jected  for  much  of  the  year. 
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5.  Under  high  May  flow  hydrology,  projected  chlorophyll 
'a'  trends  in  the  upper  channel  segments  are  similar 
to  those  under  low  flow.  In  the  shallow  segments, 
moderately  high  chlorophyll  concentrations  exist 
during  the  spring  which  then  decline  in  early  summer 
due  to  the  flushing  by  the  high  flow  in  May.  Mid¬ 
summer  peaks  of  over  100  yg/1  are  then  experienced 
in  these  shallow  segments  and  are  followed  by  sus¬ 
tained  levels  between  60  and  100  ug/1  through  the 
fall  and  early  winter. 

6.  The  patterns  of  chlorophyll  'a'  concentration  pro¬ 
jected  by  the  completely-mixed  model  are  essentially 
the  same  as  those  projected  in  the  lower  reservoir 
segments  by  the  six-segment  model. 

Lake  Model  Dissolved  Oxygen  Studies 

7.  A  biochemical  oxygen  demand/dissolved  oxygen  model 
was  developed  and  linked  to  the  six-segment  Tennessee 
Colony  Lake  eutrophication  model.  This  model,  unlike 
the  eutrophication  model,  has  not  been  calibrated 
through  testing  with  the  Lake  Livingston  eutrophi¬ 
cation  model  or  through  comparison  with  measured 
values  of  DO.  A  BOD^  rating  curve  was  developed  to 
determine  flow  related  loadings  from  the  Trinity 
River.  Daily  values  of  DO  in  each  of  the  six  segments 
were  predicted  under  both  flow  regimes. 

8.  In  order  to  evaluate  the  impact  of  the  upstream  loadings 
relative  to  the  internally-created  algal  effects, 

the  model  was  also  run  with  algal  activity  suppressed 
and  with  only  upstream  sources  of  BODg  as  inputs. 

9.  The  BOD^  rating  curve  was  adjusted  to  reflect  two 
improved  levels  of  municipal  wastewater  treatment, 
with  BODg  effluent  concentrations  of  10  mg/1  and 
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5  mg/1,  which  could  occur  in  the  future.  Daily  lake 
dissolved  oxygen  levels  were  calculated  using  the 
two  adjusted  BOD^  loading  curves,  with  and  without 
the  effect  of  algal  activity,  for  both  flow  regiir.es. 

10.  Dissolved  oxygen  levels  in  the  lake  may  be  sensitive 
to  the  degree  of  treatment  of  municipal  wastewater 
generated  in  the  Dallas/Fort  Worth  area.  Under  present 
loadings  and  low  flow  hydrology,  dissolved  oxygen 
concentrations  in  the  shallow  segments  of  the  lake 

are  projected  to  remain  close  to  or  slightly  greater 
than  saturation  throughout  the  year.  The  upper  chan¬ 
nel  segments,  however,  will  be  affected  by  the  incoming 
BOD[.  loads.  The  lake  model  indicates  that  the  June 
and  July  DO's  are  projected  to  dip  below  3  mg/1  in 
the  uppermost  channel  segments.  In  these  segments, 
recovery  in  the  late  summer  and  early  fall  is  then 
followed  by  a  less  severe  dissolved  oxygen  depression 
in  November. 

11.  Under  the  high  May  flow  hydrology  and  present  BOD,. 
loadings,  DO  values  in  June  and  July  are  projected 
to  range  from  a  low  of  2  mg/1  in  the  upper  channel 
to  a  low  of  6  mg/1  in  the  shallow  segments  in  the 
main  body  of  the  lake.  A  second  decline  to  values 
below  5  mg/1  is  projected  in  the  uppermost  channel 
sections  in  late  autumn. 

12.  Under  postulated  1979  and  long-term  BOD^  loading  con¬ 
ditions,  substantial  improvement  in  DO  levels  through¬ 
out  the  lake  is  projected.  Potential  problems  are 
only  indicated  in  the  uppermost  channel  segments, 
where  levels  slightly  below  5  mg/1  may  be  seen.  It 
should  be  kept  in  mind,  however,  that  the  model  pro¬ 
vides  only  daily  average  concentrations,  and  cannot 
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predict  the  range  of  diurnal  variation  about  that 
average.  Furthermore,  these  calculations  assume  no 
instream  nitrification.  The  results  are  to  be  con¬ 
sidered  preliminary,  as  the  BOD-DO  portion  of  the 
eutrophication  model  is  uncalibrated.  Additional 
insight  has  been  obtained  in  this  study  through  the 
use  of  the  modified  Trinity  River  model,  which  is 
discussed  below. 

Trinity  River  Model  Studies  of  Dissolved  Oxygen  in  the 

Lake  Headwaters 

13.  The  mathematical  model  of  the  Trinity  River  which 
had  been  developed  for  the  TRA  in  1974  was  recali¬ 
brated  from  its  headwaters  to  near  the  proposed  damsite 
using  comprehensive  water  quality  data  collected 
recently.  This  model  framework  is  the  basis  for  the 
evaluation  of  the  potential  impact  of  the  lake  on 

the  Trinity  River.  The  geometry  relationships  for 
the  reach  to  be  affected  by  backwater  from  the  lake 
were  modified  to  allow  the  simulation  of  dissolved 
oxygen  profiles  with  and  without  the  lake. 

14.  BOD5  and  ammonia  nitrogen  point  source  loads  in  the 
Dallas/Fort  Worth  area  were  estimated  for  1980  and 
1990  conditions.  Projected  1980  discharge  volumes 
were  associated  with  BOD^/TSS/NH-j-N  levels  of  10/15/10 
mg/1.  Projected  1990  discharge  volumes  were  asso¬ 
ciated  with  10/15/10  mg/1  and  with  5/5/3  mg/1.  Dis¬ 
solved  oxygen  profiles  in  the  upper  Trinity  River 
from  Dallas  through  the  headwaters  of  the  lake  were 
calculated  under  steady  state  low  flow,  summer  con¬ 
ditions  for  each  of  the  three  loading  conditions. 
Profiles  were  developed  with  and  without  instream 
nitrification  of  the  ammonia.  For  comparison,  profiles 
were  also  calculated  based  on  the  existing  geometry. 
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15.  Projected  1980  and  1990  DO  profiles  in  the  modified 
Trinity  River  model  indicate  a  reduction  of  between 
1.5  and  3.0  mg/1  in  the  headwaters  of  the  proposed 
lake  if  the  dam  were  constructed  as  compared  with 
conditions  in  the  river  if  the  dam  were  not  constructed. 
These  results  are  projected  to  occur  during  low  flow 
summer  conditions. 

16.  Whether  nitrification  will  occur  in  the  Trinity  River 
downstream  of  advanced  waste  treatment  plants  is  uncer¬ 
tain.  The  results  indicate  a  signficant  difference 

in  the  calculated  DO  in  the  lake  headwaters  with  and 
without  instream  nitrification.  The  1990  loadings 
at  effluent  concentrations  of  10  mg/1  BOD^  and  10  mg/1 
NH^-N  would  cause  a  DO  depression  to  less  than  1  mg/1 
in  the  headwaters  with  instream  nitrification  and 

4  mg/1  without.  Further,  for  advanced  treatment  to 

5  mg/1  BOD5  and  3  mg/1  NH^-N ,  these  DO  levels  would 
be  approximately  3  mg/1  with  instream  nitrification 
and  4.5  mg/1  without. 

Trinity  River  Model  Studies  of  Dissolved  Oxygen  Downstream 

of  the  Dam 

17.  An  evaluation  of  dissolved  oxygen  conditions  down¬ 
stream  of  the  proposed  dam  was  prepared.  A  one-hundred 
mile  reach  of  the  river  was  examined  using  a  section 

of  the  previously  developed  TRA  model.  In  order  to 
investigate  the  general  recovery  characteristics  of 
the  river  to  loads  carried  in  the  dam  release  water, 
six  combinations  of  dissolved  oxygen  (DO)  and  ultimate 
oxygen  demand  (UOD)  concentrations  in  the  release 
water  were  applied  under  each  of  the  two  flow  con¬ 
ditions.  This  provided  an  indication  of  the  sensi¬ 
tivity  of  the  river  DO  profile  to  the  DO  and  UOD  in 
the  release  water. 


18.  The  Trinity  River  downstream  of  the  damsite  appears 
to  recover  slowly  from  dissolved  oxygen  depression. 
Therefore,  the  river  will  be  sensitive  to  the  amount 
of  oxygen-demanding  material  present  in  the  release 
water  and  to  the  dissolved  oxygen  deficit  of  the 
water.  If  the  oxygen  demand  in  the  water  is  low, 
such  as  an  ultimate  oxygen  demand  (UOD)  of  5  mg/1, 
the  minimum  DO  in  the  river  will  be  at  the  release 
point  itself.  In  this  case,  to  meet  a  stream  standard 
of  5  mg/1,  the  release  water  DO  concentration  should 
be  5  mg/1  or  slightly  greater.  If  higher  levels  of 
oxygen  demanding  material  are  present,  such  as  a  UOD 
of  15  mg/1  or  greater,  release  water  DO  levels  in 
excess  of  5  mg/1  would  be  needed  to  protect  against 
violation  of  the  stream  standard. 
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CHAPTER  II 
INTRODUCTION 

This  report  summarizes  Phase  II  of  a  water  quality 
study  of  Lake  Livingston  and  the  proposed  Tennessee  Colony 
Lake.  The  sites  of  these  reservoirs  are  located  on  the 
main  stem  of  the  Trinity  River,  Texas.  The  previous  Phase  I 
e-fort^  investigated  the  projected  eutrophication  condi¬ 
tions  of  the  final  configuration  of  the  Tennessee  Colony 
Lake  with  a  conservation  pool  elevation  of  275  feet.  The 
present  Phase  II  effort  has  projected  conditions  for  an 
interim  pool  configuration  with  a  conservation  pool  eleva¬ 
tion  of  240  feet.  Plate  1  presents  a  map  of  this 
configuration. 

The  following  sections  of  this  chapter  will  provide 
introductory  material  briefly  reviewing  the  Phase  I  work 
and  outlining  the  purpose  and  scope  of  the  Phase  II  work. 

A.  Summary  of  Phase  I  Effort. 

In  the  initial  effort,  a  eutrophication  model  was 
developed  for  Lake  Livingston,  an  existing  reservoir  approxi¬ 
mately  213  river  miles  downstream  of  the  location  of  the 
proposed  Tennessee  Colony  Lake  dam.  Chlorophyll  'a'  was 
selected  as  an  indicator  of  algal  growth.  This  model, 
which  tracked  suspended  solids  and  nutrients  through  the 
lake  during  an  entire  year  in  order  to  calculate  chloro¬ 
phyll  'a'  levels,  was  calibrated  using  1975  data  and  vali¬ 
dated  using  1976  data.  This  procedure  produced  appropriate 
values  for  various  kinetic  parameters  such  as  decay  rates, 
settling  rates  and  dispersion  rates.  A  similar  model 
structure  was  then  developed  for  the  final  configuration 
of  the  proposed  Tennessee  Colony  Lake  in  which  the  Lake 
Livingston  kinetic  parameters  were  used.  Yearly  profiles 
of  chlorophyll  'a'  were  then  projected  for  Tennessee  Colony 
Lake  for  two  distinct  flow  regimes. 
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MAP  OF  INTERIM  POOL  CONFIGURATION 
PROPOSED  TENNESSEE  COLONY  LAKE 


B .  Purpose  of  Phase  II  Effort. 

In  the  Phase  II  study,  a  similar  model  was  constructed 
for  a  smaller  interim  pool  for  which  the  kinetic  parameters 
developed  in  Phase  I  were  used.  This  model  was  used  to 
project  water  quality  conditions  in  the  interim  pool. 

As  in  the  previous  effort,  suspended  solids  and  nutrients 
were  tracked  through  the  reservoir  and  concentrations  of 
chlorophyll  'a'  throughout  the  year  were  calculated.  In 
addition,  calculations  of  dissolved  oxygen  levels  through¬ 
out  the  year  were  made.  This  portion  of  the  model  is 
uncalibrated  as  it  was  not  included  in  the  original  Lake 
Livingston  model. 

A  second  area  of  concern  in  the  Phase  II  effort  was 
the  potential  impact  of  the  Tennessee  Colony  Lake  on  water 
quality  in  the  Trinity  River.  Two  river  dissolved  oxygen 
models  were  developed  to  investigate  this  impact:  one 
for  the  upper  Trinity,  including  the  headwaters  of  the 
Tennessee  Colony  Lake,  and  one  for  the  portion  of  the  river 
extending  100  miles  downstream  of  the  dam.  These  models 
were  applied  primarily  to  investigate  projected  dissolved 
oxygen  conditions  for  conditions  for  steady  state,  low 
flow  summer  conditions.  This  low  flow  summer  period 
generally  represents  the  most  critical  period  of  extended 
depressed  dissolved  oxygen  due  to  higher  reaction  rates, 
reduced  DO  saturation,  less  dilution  water,  and  lowered 
stream  velocities. 

Chapter  III  of  this  report  will  discuss  the  results 
of  the  Tennessee  Colony  eutrophication  model;  Chapter  IV 
will  discuss  the  results  of  the  two  Trinity  River  water 
quality  models. 
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CHAPTER  III 

TENNESSEE  COLONY  LAKE  EUTROPHICATION 


CHAPTER  III 

TENNESSEE  COLONY  LAKE  EUTROPHICATION 


A  lake  eutrophication  model  was  constructed  for  the 
interim  configuration  of  the  proposed  Tennessee  Colony 
Lake.  As  in  the  models  developed  for  Lake  Livingston  and 
the  final  configuration  of  Tennessee  Colony  Lake,  which 
are  discussed  in  the  Phase  I  report, ^  chlorophyll  'a' 
is  used  as  an  indicator  of  algal  growth,  a  chief  measure 
of  eutrophication.  Projections  of  chlorophyll  levels  throug 
out  the  year  are  presented  for  two  annual  flow  regimes. 

In  addition,  simulations  of  BOD-DO  conditions  in  the 
interim  pool  were  performed.  This  additional  BOD-DO  model 
was  constructed  only  for  Tennessee  Colony  Lake  and  was 
not  included  in  the  previous  1975  and  1976  Lake  Livingston 
calibration  and  validation. 

The  following  sections  of  this  chapter  will  describe 
the  computational  framework  of  the  models,  review  the  geo¬ 
metry,  segmentation  and  transport  mechanisms  used  for  pro¬ 
jection  purposes,  and  present  the  results  of  projections 
made  using  the  models. 

A.  Description  of  Model  Framework. 

The  eutrophication  model  employed  in  this  study  was 
developed  as  a  practical  engineering  tool  to  allow  for 
the  inclusion  of  the  major  conditions  influencing  phyto¬ 
plankton  growth.  The  model  utilizes  a  growth  relation¬ 
ship  which  includes  the  major  physiological  factors  which 
in  most  cases  dominate  the  algal  growth  relationship. 

Growth  is  related  to  the  temperature  of  the  water  body, 
the  available  sunlight,  suspended  solids  concentrations, 
and  the  concentrations  of  the  major  nutrient  forms  which 
can  be  directly  utilized  by  algae.  Each  of  these  control¬ 
ling  elements  interacts  in  defining  the  algal  growth  rate. 
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Nutrient  concentrations  are  in  turn  related  to  changes 
in  incoming  loads  and  to  algal  uptake. 

The  eutrophication  model  is  simply  a  computational 
framework  within  which  these  complex  interrelationships 
can  be  handled  simultaneously.  The  basis  of  the  model 
is  a  series  of  conservation  of  mass  equations  which  relate 
the  variables  of  concern  to  each  other  and  also  to  the 
boundary  loading  conditions.  The  mass  conservation  equa¬ 
tions  which  result  from  this  formulation  are  partial  dif¬ 
ferential  equations  in  the  time  and  space  variables.  In 
order  to  implement  the  solution  of  such  equations  on  a 
computer,  it  is  convenient  to  express  them  in  terms  of 
finite  differences.  For  the  spatial  variables,  this  corre¬ 
sponds  to  dividing  the  water  body  into  a  series  of  segments 
or  cells  which  are  chosen  so  that  the  assumption  of  spatial 
homogeneity  within  each  segment  is  reasonable.  A  typical 
conservation  of  mass  equation  used  in  the  model  for  con¬ 
centration  c^  of  substance  i  in  segment  j  has  the  general 
form: 


dc  •  • 


dt 
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kj 


(cik  ‘ 
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where  Vj  is  the  segment  volume;  is  the  net  source  of 

substance  i  in  segment  j;  is  the  bulk  rate  of  transport 

of  c . .  into  and  c. .  out  of  segment  j  for  all  segments  k 

IK  1  J 

adjacent  to  segment  j;  and  Q^j  is  the  net  advective  flow 
rate  between  segments  k  and  j.  Numerical  integration  of 
these  equations  gives  the  daily  distribution  of  the  concen¬ 
trations  in  each  of  the  spatial  segments  of  the  model. 

A  detailed  description  of  model  kinetics  has  been 
provided  as  Appendix  B  to  the  Phase  I  report. 
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B.  Model  Geometry. 


The  Tennessee  Colony  Lake  model  for  the  interim  pool 
is  shown  in  Plate  1  and  schematically  in  Figure  1.  The 
model  includes  the  area  of  the  proposed  impoundment  that 
would  be  covered  at  the  conservation  pool  elevation  of 
240  feet  above  sea  level,  some  15,100  acres.  The  pool 
will  contain  a  volume  of  109,900  acre  feet  or  4,787  million 
cubic  feet.  Lake  level  is  assumed  to  remain  constant 
throughout  the  year;  that  is,  flow  out  equals  flow  in. 

Two  model  configurations  have  been  developed  to  repre¬ 
sent  two  alternative  approaches  to  the  transport  and  mixing 
characteristics  of  the  proposed  impoundment.  The  first 
configuration  consists  of  six  segments.  Three  segments 
correspond  to  the  deep,  existing  river  channel  areas  of 
the  impoundment;  three  segments  are  used  to  describe  the 
large,  shallow  areas  adjacent  to  the  deeper  river  channel. 
The  second  configuration  assumes  the  impoundment  to  be 
a  single  completely-mixed  volume,  with  total  lake  volume 
and  average  lake  depth. 

The  geometries  of  the  model  segments  for  both  configu¬ 
rations  are  presented  in  Table  1.  Segment  interconnections 
are  defined  by  characteristic  lengths,  dispersion  coeffi¬ 
cients,  and  interfacial  cross-sectional  areas,  also  pre¬ 
sented  in  Table  1. 

C.  Transport  Regimes. 

A  basic  requirement  for  a  modeling  effort  based  on 
conservation  of  mass  is  an  adequate  representation  of  the 
mass  transport  mechanisms  in  the  study  area.  The  represen¬ 
tation  chosen  is  based  on  the  advection-dispersion  for¬ 
mulation  of  mass  transport.  Advective  mass  transport  is 
accomplished  by  the  average  unidirectional  net  motions 
of  the  water.  Dispersive  mass  transport  is  accomplished 
by  the  mixing  motions  of  the  water  body  such  as  the  smaller 
scale  circulations. 


FIGURE  I 

TENNESSEE  COLONY  LAKE  INTERIM  POOL  CONFIGURATION 

MODEL  SCHEMATIC 


TABLE  3 


TENNESSEE  COLONY  LAKE  INTERIM  POOL  GEOMETRY 


SEGMENT  GEOMETRY 


SEGMENT 

VOLUME 

(106ft3) 

SURFACE 

AREA 

(103ft2) 

MEAN 

DEPTH 

(ft) 

1 

491.9 

18.0 

27.3 

2 

40.7 

2.22 

18.3 

3 

618.1 

42.5 

14.5 

Channel  Segment 
Subtotal 

1,150.7 

62.7 

(16.4) 

4 

807.4 

115.4 

7.0* 

5 

2,124.1 

303.4 

7.0* 

6 

705.0 

176.2 

4.0* 

Shallow  Segment 
Subtotal 

3,636.5 

595.0 

(6.1) 

Reservoir  Total 

4,787.2 

657.7 

(7.3) 

INTERFACE  GEOMETRY 

DEPTH  AT 


WIDTH  INTERFACE 
INTERFACE  (ft)  (ft) 


180  42.0 

1- 2  277  23.5 

2- 3  280  18.7 

3- 0  112  9.9 

41,000  12.5 

1- 5  48,000  12.5 

2- 6  8,000  7.0 


CROSS- 

SECTIONAL 

AREA 

(ft2) 

LENGTH 

I 

(ft) 

LENGTH 

J 

(ft) 

7,564 

48,000 

48,000 

6,502 

48,000 

12,000 

5,232 

12,000 

147,500 

1,109 

147,500 

147,500 

512,500 

229** 

3,000 

600,000 

229** 

9,000 

56,000 

279*** 

11,400 

*  Estimated 

♦*  Average  of  0-1  and  1-2  width 
•♦♦Average  of  1-2  and  2-3  width 


Flow  regimes  into  Tennessee  Colony  Lake  were  deter¬ 
mined  based  on  analyses  of  flow  estimates  prepared  by  the 
Corps  of  Engineers'  The  Corps  of  Engineers  has  presented 
forty-seven  years  of  monthly  flows  at  Tennessee  Colony 
adjusted  for  existing  and  authorized  projects  in  the  Trinity 
River  basin.  Several  recurring  patterns  are  evident  upon 
examination  of  the  annual  hydrologic  record.  First,  the 
summer  period  of  July,  August,  and  September  typically 
has  low  flows,  generally  less  than  650  cfs  at  Tennessee 
Colony.  Second,  annual  hydrology  can  be  described  on  the 
basis  of  the  magnitude  of  May  flows.  The  characterization 
of  an  annual  record  as  having  either  a  high  or  a  low  May 
flow  is  a  prominent  feature  of  the  data.  Figures  2  and  3 
present  examples  of  years  which  have  high  monthly  average 
flows  in  May  and  those  years  having  low  flows  in  May. 

Figure  2  presents  Trinity  River  flows  at  Tennessee  Colony 
for  the  six  years  in  the  record  with  the  lowest  average 
monthly  May  flows.  Figure  3  presents  the  Trinity  River 
flows  at  Tennessee  Colony  for  the  six  years  in  the  record 
with  the  highest  May  flows.  In  each  instance  logarithmic 
averages  of  the  monthly  flows  are  calculated  to  produce 
a  composite  annual  profile  which  is  indicated  on  the  figures. 
The  log  averaging  process  is  used  since  the  intent  of  the 
analysis  is  to  show  patterns  rather  than  an  average  flow. 

The  log  average  tends  to  damp  out  extremes  in  the  data 
and  thus  is  more  indicative  of  average  trends. 

These  two  flow  patterns,  that  is,  years  with  high 
flows  in  May  and  years  with  low  flows  in  May,  are  repre¬ 
sentative  of  the  majority  of  annual  patterns  in  the  record 
and  present  limits  which  are  used  to  bracket  flows  which 
can  reasonably  be  expected  in  the  future.  The  log  averages 


highest  mat  flow 


F,GURE  3  ETADC 

MONTHLY  AVERAGE  FLOWS  AT  TENNESSEE  COLONY,  HIGH  MAY  FLOW  YEARS 


of  the  monthly  flows  are  used  for  projection  purposes  in 
the  Tennessee  Colony  Lake  model,  as  discussed  in  subse¬ 
quent  sections  of  this  chapter. 

Flow  is  routed  through  the  main  channel  segments  of 
the  six-segment  model,  and  in  and  out  of  the  completely- 
mixed  model.  Constant  segment  and  reservoir  volumes  are 
assumed,  i.e.,  whatever  flow  volume  enters  from  up  river 
is  discharged  at  the  dam. 

One  measure  of  sensitivity  to  flow  conditions  is 
hydraulic  detention  time.  Hydraulic  detention  time,  defined 
as  volume  divided  by  flow,  is  the  average  time  required 
for  water  to  pass  through  a  volume.  For  conservative 
substances,  approximately  two  detention  times  are  necessary 
for  a  completely-mixed  volume  to  approach  a  new  equilibrium 
condition.  Hydraulic  detention  times  based  upon  the  total 
interim  pool  volume  of  Tennessee  Colony  Lake  range  from 
over  nine  months  at  low  flows  below  200  cfs  to  less  than 
six  days  at  flows  above  10,000  cfs.  Average  low  flow  year 
detention  time  is  approximately  2.5  months;  high  May  flow 
year  detention  time  is  on  the  order  of  13  days.  At  peak 
high  flow  during  the  high  May  hydrology  of  28,300  cfs, 
hydraulic  detention  time  decreases  to  approximately  two 
days.  These  numbers  can  be  understood  as  follows;  at 
low  flows,  even  large  changes  in  influent  concentrations 
will  have  minor  effects  on  the  concentration  of  the  com¬ 
pletely-mixed  volume  for  most  time  frames;  at  higher  flows, 
changes  in  influent  concentrations  have  increasingly  greater 
impacts  on  completely-mixed  volume  concentrations. 

Horizontal  dispersion  coefficients  are  used  to  provide 
mixing  between  the  deep  channel  segments  and  the  shallow 
segments  of  the  flats.  A  range  of  transport  coefficients 
was  examined  since  the  Tennessee  Colony  Lake  interim  period 
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pool  configuration  and  volume  are  quite  different  from 
those  of  the  Lake  Livingston  reservoir  which  was  used  for 
model  calibration.  Results  were  produced  using  a  value 
of  0.1  mi  /day,  the  value  used  in  the  Lake  Livingston  model 
and  the  Phase  I  Tennessee  Colony  model.  The  sensitivity 
of  the  model  results  to  this  value  was  examined  by  repro¬ 
ducing  calculated  concentrations  using  a  value  of  0.002 
mi  /day.  This  value  is  consistent  with  those  typical  of 
small  and  medium-sized  rivers.  Comparison  of  the  results 
obtained  using  this  range  of  values,  plus  the  results  from 
the  completely-mixed  volume  assumption,  should  provide 
an  understanding  of  the  sensitivity  of  model  projections 
to  transport  coefficients. 

D.  Loading  Conditions. 

For  the  purposes  of  this  study,  all  nutrient,  BOD 
and  suspended  solids  loadings  to  Tennessee  Colony  Lake 
are  assumed  to  be  associated  with  Trinity  River  inflows 
to  the  impoundment.  Other  sources,  such  as  Richland  Creek, 
Chambers  Creeks,  and  Tehuacana  Creek,  join  the  Trinity 
River  above  the  main  body  of  the  interim  pool.  No  direct 
point  or  nonpoint  source  discharges  are  postulated  at  the 
present  time. 

An  in-depth  analysis  of  constituent  concentrations 
in  the  Trinity  River  at  Rosser,  approximately  60  miles 
upstream  from  the  impoundment,  was  performed  as  part  of 
the  Phase  I  effort.  The  relationships  between  flow  and 
concentration  were  examined  for  chlorophyll  'a',  suspended 
solids,  nitrogen  forms,  phosphorus  forms,  and  BOD^;  with 
the  exception  of  chlorophyll  'a',  rating  curves  were 
developed  for  all  constituents  which  serve  to  provide  values 
of  concentration  for  given  flows.  Daily  loadings  via  the 


23 


Trinity  River  inflow  were  input  to  the  Tennessee  Colony 
Lake  model  by  selecting  the  appropriate  concentrations 
from  these  curves. 

1.  Rating  Curves.  Analysis  of  the  concentration 
data  indicated  a  strong  positive  correlation  between 
suspended  solids  and  flow  and  strong  negative  correlation 
between  BODj,  nitrogen  series,  and  phosphorus  series  con¬ 
centrations  and  flow.  Little  variation  in  chlorophyll 
'a'  levels  was  shown.  The  trends  of  the  chemical  consti¬ 
tuents  are  shown  in  Figures  4  and  5. 

Figure  4  presents  values  of  total  nitrogen,  ammonia 
nitrogen,  nitrate  plus  nitrite  nitrogen,  and  organic  nitro 
gen  versus  river  flow  at  Rosser.  Figure  5  presents  values 
of  total  phosphorus,  BOD^,  and  total  suspended  solids  ver¬ 
sus  flow  for  the  same  location.  Data  are  from  published 
U.S.  Geological  Survey  (USGS)  records  for  the  period  1971 
through  1975  ^.  Where  lines  are  indicated  on  the  figures 
the  lines  are  trend  lines,  not  regressions. 

Total  nitrogen  and  total  phosphorus  concentrations 
exhibit  similar  behavior  with  flow.  High  total  nutrient 
concentrations  are  characteristic  of  low  river  flows: 
concentrations  tend  to  decrease  as  flow  increases, 
apparently  due,  at  least  in  part,  to  dilutional  effects. 
Further,  there  appears  to  be  a  possible  background  level 
concentration  even  at  higher  flows. 

These  patterns  can  be  explained  as  follows:  at  low 
flow,  the  major  portion  of  Trinity  River  flow  is  composed 
of  treated  effluent  from  the  Dallas-Fort  Worth  Metroplex; 
as  flow  increases  due  to  land  runoff,  concentrations 
decrease  due  to  dilution  with  lower  concentration  waters, 
until  at  high  flows,  the  river  concentrations  approach 
values  typical  of  those  expected  in  nonpoint  source  runoff 
for  east  Texas (*'5,6) ^  In  particular,  at  Rosser,  total 
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FIGURE  4 

MEASURED  NITROGEN  SERIES  CONCENTRATIONS  VERSUS  RIVER  FLOW 

AT  ROSSER,  TEXAS 


FIGURE  5 

TOTAL  PHOSPHORUS,  BODb,  AND  TOTAL  SUSPENDED  SOLIDS 
CONCENTRATIONS  VERSUS  RIVER  FLOW  AT  ROSSER,  TEXAS 


nitrogen  concentrations  range  from  a  high  value  of  approxi¬ 
mately  15  mg/1  at  low  flows  to  concentrations  of  4  to  2 
mg/1  at  moderate  flows  (2,000  -  6,000  cfs)  ,  to  approxi¬ 
mately  1.2  mg/1  at  flows  above  10,000  cfs.  Total  phosphorus 
concentrations  range  from  a  high  of  about  9  mg/1  at  low 
flows,  to  concentrations  of  4  to  0.7  mg/1  at  moderate  flows, 
to  approximately  0.5  mg/1  at  flows  above  10,000  cfs. 

The  individual  nitrogen  species  are  also  of  interest 
both  for  an  understanding  of  what  is  occurring  in  the 
system  and  because  a  breakdown  by  nitrogen  species  is 
necessary  for  model  input.  At  Rosser,  ammonia  nitrogen 
contributes  the  major  fraction  of  the  total  nitrogen  con¬ 
centration,  and  consequently,  the  ammonia  nitrogen  and 
the  total  nitrogen  concentrations  range  from  about  13  mg  N/l 
at  low  flows  to  concentrations  of  2.0  to  0.5  mg  N/l  at 
moderate  flows  to  approximately  0 . 3  mg  N/l  at  higher  flows. 

Conversely,  nitrate  plus  nitrite  nitrogen  and  organic 
nitrogen  concentrations  exhibit  no  such  flow  dependence 
at  Rosser.  Figure  4  presents  nitrate  plus  nitrite  nitrogen 
at  Rosser  as  a  function  of  flow  for  those  data  having 
synoptic  measurements  for  ammonia  and  organic  nitrogen. 

The  data  are  fairly  constant  at  all  flows,  and  have  a  mean 
value  of  0.8  mg  N/l.  Figure  4  also  presents  total  organic 
nitrogen  at  Rosser  as  a  function  of  flow.  While  more 
scattered  than  the  nitrate  plus  nitrite  data,  the  shape 
is  similar,  with  no  discernible  flow  related  effects 
present.  The  mean  value  for  total  organic  nitrogen  at 
Rosser  is  0.9  mg  N/l. 

Figure  5  presents  BOD^  measurements  at  Rosser  versus 
river  flow.  While  data  scattering  is  evident,  a  statis¬ 
tical  grouping  of  measurements  indicates  a  relatively  con¬ 
stant  mean  BOD&  concentration  of  12.8  mg/1  until  river 
flows  exceed  3,000  cfs.  Mean  BOD^  concentrations  decrease 
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markedly  above  3,000  cfs.  BOD^  data  at  flows  above  5,000 
cfs  are  sparse. 

It  should  be  recognized  that  these  BOD^  data  relate 
to  the  total  oxygen  demand  of  a  sample  during  a  5-day 
incubation  period.  As  such,  the  measurements  include  the 
oxygen  utilized  by  any  of  several  mechanisms,  including 
bacterial  decomposition  and  algal  respiration.  As  bacterial 
decomposition  of  carbonaceous  material  is  the  primary  process 
of  interest  in  the  BOD^  measurement,  the  effects  of  algal 
respiration  should  be  discounted.  When  measurements  of 
algal  biomass  are  made  concurrent  with  BOD^  measurements, 
it  is  possible  to  estimate  the  fraction  of  the  measured 
BOD5  which  is  due  to  algal  respiration.  The  measured  BOD_ 
values  can  then  be  adjusted  to  reflect  only  bacterial 
activity.  No  such  concurrent  measurements  are  available 
in  the  present  study;  these  BOD5  values  therefore  are 
not  totally  reflective  of  the  available  BODj.  entering  the 
Tennessee  Colony  Lake. 

Figure  5  presents  total  suspended  solids  concentration 
versus  river  flow  at  Rosser.  These  data  are  daily  measure¬ 
ments  covering  the  period  1971  to  1975  and  are  available  from 
the  Texas  Natural  Resources  Information  System.  The  data 
are  variable,  but  show  increasing  concentration  with  flow 
up  to  some  point  above  which  concentration  appears  to 
decrease . 

2.  Boundary  Conditions.  The  nutrient,  BOD^,  and 
total  suspended  solids  rating  curves  are  applied  to  the 
high  or  low  May  annual  hydrologic  regimes  to  provide  loading 
information  to  the  eutrophication  model.  This  is  done 
via  model  boundary  conditions.  For  a  particular  daily 
flow,  an  appropriate  concentration  is  selected  from  the 
rating  curve.  This  concentration  is  applied  at  the  model 
boundary  to  input  the  proper  mass  of  material  to  the  model. 
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An  exception  to  this  procedure  is  the  determination 
of  chlorophyll  'a'  boundary  conditions.  Inflow  levels 
were  set  to  a  constant  5  pg/1  for  all  flow  conditions. 

It  was  determined  during  the  Lake  Livingston  reservoir 
calibration  and  validation  procedure  that  data  regarding 
higher  chlorophyll  'a'  concentrations  in  the  Lake  Livingston 
inflow  notwithstanding,  a  constant  boundary  condition  would 
be  adequate  to  reproduce  observed  reservoir  chlorophyll 
'a'  concentrations.  This  results  from  a  combination  of 
effects.  First,  because  of  the  magnitude  of  the  algal 
growth  and  death  rates,  phytoplankton  growth  kinetics  in 
the  reservoir,  rather  than  chlorophyll  'a'  influent  con¬ 
centrations,  determine  algal  levels  in  the  lake  for  most 
flow  conditions.  Second,  with  regard  to  a  strict  mass 
balance  of  materials,  the  mass  of  chlorophyll  'a'  entering 
the  reservoir  is  generally  small  and  does  not  appreciably 
alter  reservoir  concentrations  due  to  the  large  volume 
of  the  lake.  This  relative  insensitivity  to  the  chloro¬ 
phyll  'a'  boundary  condition  is  fortuitous  since  it  ob¬ 
viates  the  need  for  precise  chlorophyll  'a'  boundary 
concentration  determination. 

E.  Initial  Conditions. 

It  is  necessary  in  a  time  variable  simulation  to  supply 
initial  concentration  values  for  each  system  variable  in 
each  model  segment  as  a  starting  point  for  computations. 

The  choice  of  these  initial  conditions  will  have  an  effect 
on  subsequent  computations  until  such  time  as  the  initial 
conditions  have  "washed  out"  of  the  system.  These  are 
called  transient  effects.  In  terms  of  the  previous  dis¬ 
cussion  of  hydraulic  detention  time,  approximately  two 
detention  times  are  required  for  transient  effects  to 
be  minimized. 


Ideally,  initial  conditions  are  either  chosen  from 
observed  data  or  the  model  is  run  for  several  years  of 
simulation  to  wash  out  assigned  initial  conditions.  Since 
in  the  present  case,  the  latter  method  requires  additional 
assumptions  as  to  the  repetitiveness  of  annual  hydrologies, 
model  initial  conditions  in  all  cases  were  assigned,  as 
best  estimates,  from  the  rating  curves. 

F .  Model  Projections. 

All  of  the  above  factors  were  incorporated  into  the 

eutrophication  model  of  the  interim  pool  of  the  proposed 

Tennessee  Colony  Lake.  The  model  was  used  to  project  annual 

water  quality  trends  under  the  two  idealized  flow  regimes. 

The  single-segment,  completely-mixed  model  projected  values 

of  chlorophyll  'a',  an  indicator  of  algal  biomass,  under 

the  high  May  and  low  May  flow  regimes.  The  six-segment 

model  projected  values  of  chlorophyll  'a'  under  both  flow 

regimes,  once  using  a  moderate  horizontal  dispersion  rate 

(0.1  mi  /day)  and  then  using  a  lower  bound  value 
,  2 

(0.002  mi  /day).  In  addition,  projections  of  dissolved 
oxygen  were  produced  using  the  six-segment  configuration 
under  various  flow  and  loading  conditions.  These  projec¬ 
tions  examine  the  expected  impoundment  responses  for  changes 
in  hydrologic  regimes  and  the  sensitivity  of  the  model 
to  assumptions  regarding  transport  within  the  lake.  It 
must  be  emphasized  that  dissolved  oxygen  projections  from 
the  eutrophication  model  are  only  preliminary  estimates 
since  dissolved  oxygen  model  kinetics  are  uncalibrated 
for  the  lake  model.  The  estimates  can  however  be  useful 
in  the  evaluation  of  the  proposed  project  and  in  the  analysi 
of  specific  designs  and  operation  procedures.  Table  2 
summarizes  the  model  projections  produced  for  the  interim 
pool. 
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SUMMARY  OF  INTERIM  TENNESSEE  COLONY  LAKE  MODEL  RUNS 
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Note:  No  nitrification  of  ammonia  nitrogen. 

*  -Upstream  loads  only-  indicates  suppression  of  algal  impact  on  production  and  demand  of  oxygen 


1.  Chlorophyll  1  a 1 .  Chlorophyll  'a',  as  a  measure 
of  phytoplankton  biomass,  is  the  model  variable  by  which 
comparisons  of  lake  trophic  state  are  made.  It  is  a  direct 
measure  of  the  quantity  of  photosynthetic  plant  pigments 
present,  and  can  be  related  directly  to  oxygen  production. 
Discussions  of  chlorophyll  'a'  as  an  indicator  of  trophic 
level,  water  quality,  and  beneficial  water  uses  have  been 
presented  in  previous  reports^'^.  Many  studies  have 
used  chlorophyll  ’a’  as  an  indicator  of  primary  productivi  ty 
or  lake  eutrophic  state ^ ' 10 ' 11 ' 12 ' .  Interpretation 
of  the  impact  of  algal  levels  is  somewhat  subjective. 

There  is  no  universal  standard  for  acceptable  levels  of 
chlorophyll  ’a1  and  the  impact  of  this  parameter  must  be 
tied  to  local  water  use  objectives.  It  remains  for  the 
appropriate  planning  and  regulatory  agencies  to  determine 
desirable  concentration  levels. 


As  indicated  in  Table  2,  six  projections  of  chloro¬ 
phyll  'a'  were  performed.  Two  projections  are  based  on 
the  completely-mixed  configuration  and  the  remaining  four 
for  the  six-segment  model. 

a.  Completely-Mixed  Configuration.  Figure  6 
presents  chlorophyll  'a'  concentration  projections  for 
the  interim  pool  configuration  of  Tennessee  Colony  Lake 
assuming  the  entire  impoundment  is  completely  mixed. 
Calculations  are  presented  for  both  low  and  high  flow 
hydrologies . 

Chlorophyll  'a'  concentration  levels  between  40  to 
80  yg/1  are  projected  for  most  of  the  low  flow  year. 
Somewhat  different  algal  behavior  is  projected  for  the 
high  flow  year  hydrology.  Impoundment  chlorophyll  'a' 
concentrations  are  projected  to  increase  early  in  the  year 
and  then  decrease  during  the  high  flow  period.  The  lack 
of  algal  growth  during  the  high  flow  portion  of  the  year 
is  due  to  a  combination  of  light  limitation  and  detention 
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FIGURE  6 

CHLOROPHYLL  'a'  CONCENTRATION  PROJECTIONS 
COMPLETELY  MIXED  CONFIGURATION  -  EFFECT  OF  FLOW  REGIME 


time  effects.  Flow  related  increases  in  influent  suspended 
sediments  cause  an  increase  in  impoundment  suspended  solids 
concentrations,  causing  a  corresponding  decrease  in  available 
light  for  phytoplankton  growth.  This  decrease  causes  a 
decline  in  the  phytoplankton  growth  rate,  which  of  itself 
or  in  combination  with  the  flow  through  the  impoundment 
removing  material  from  the  impoundment,  prevents  net  in¬ 
creases  in  algal  levels.  This  occurs  even  though  nutrients 
are  plentiful. 

As  flow  decreases,  around  day  160,  light  penetration 
increases,  increasing  the  growth  rate;  impoundment  residence 
time  also  increases  due  to  decreasing  flow,  thus  allowing 
time  for  growth  to  occur.  The  result  is  a  rapid  increase 
in  chlorophyll  'a'  concentrations  to  levels  over  150  ug/1. 

At  this  point,  available  inorganic  nutrients  are  exhausted 
and  the  algal  population  declines  to  a  level  that  can  be 
sustained  by  nutrient  inflow  to  the  impoundment  and  by 
nutrient  recycle.  Nutrient  recycle  involves  the  respiration 
and  death  of  algal  cells  to  non-living  organic  nutrients, 
and  the  subsequent  conversion  to  inorganic  nutrients  avail¬ 
able  for  further  algal  growth.  This  combination  of  inflow 
and  recycle  is  calculated  to  allow  population  levels  of 
40  to  60  yg/1  chlorophyll  'a'  to  be  sustained  during  low 
flow  periods. 

b.  Multi-Segment  Configuration.  In  an  effort 
to  address  the  sensitivity  of  model  calculations  to  assump¬ 
tions  regarding  impoundment  dispersional  transport,  the 
multi-segment  model  configuration  is  used  for  projections 
and  three  deep  areas.  The  deep  areas  correspond  to  the 
existing  river  channel  running  through  the  proposed  lake 
area.  River  flow  is  assumed  to  be  confined  to  these  seg¬ 
ments.  Transport  between  deep  and  shallow  segments  is 
limited  to  dispersional  mixing. 
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Two  values  of  the  horizontal  dispersion  coefficient 
are  considered  in  the  analysis.  The  first  calculation 
employs  a  value  consistent  with  that  used  in  the  Lake 
Livingston  model  calibration  and  validation,  and  for  the 
projections  performed  for  the  original  Tennessee  Colony 
Lake  study.  The  second  calculation  reduces  the  dispersion 
coefficient  to  a  level  associated  with  small  or  medium-sized 
rivers.  The  results  of  these  computations  are  shown  in 
Figure  7  for  the  low  May  flow  hydrology  and  in  Figure  8 
for  the  high  May  flow  hydrology. 

For  low  flow  hydrology,  chlorophyll  'a'  concentrations 
are  projected  to  range  from  5  to  100  ug/1.  Considering 
the  higher  dispersion  case,  concentrations  in  Segments  2 
and  3,  the  upper  and  middle  main  stem  segments,  remain 
relatively  low  until  midsummer  (day  150  to  180) ,  at  which 
time  substantial  growth  is  projected  to  occur.  This  be¬ 
havior  is  due  to  a  combination  of  flow  related  effects: 
light  limitation  of  the  algal  growth  rate  due  to  flow 
related  suspended  solids,  and  hydraulic  detention  time  of 
the  segment.  No  net  concentration  increase  is  projected 
in  Segment  3  until  such  time  as  the  flow  through  the  segment 
is  reduced  to  very  low  levels. 

The  situation  in  the  rest  of  the  impoundment  is  quite 
different.  Rapid  growth  is  projected  in  the  remainder  of 
the  reservoir  with  chlorophyll  'a*  levels  remaining  between 
20  and  100  pg/1.  In  general,  growth  occurs  in  the  shallow 
areas.  Segments  4,  5  and  6,  and  the  algae  are  dispersed 
to  the  deeper  channel  segments.  The  shallow  areas  utilize 
the  nutrients  which  are  brought  into  the  impoundment  with 
the  flow  and  moved  to  the  shallows  by  dispersion.  The  effect 
of  the  midsummer  bloom  in  Segments  2  and  3  is  quite  clear 
in  the  remainder  of  the  reservoir.  The  bloom  converts  inor¬ 
ganic  nitrogen  to  algal  cells  in  Segments  2  and  3,  leaving 
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CHLOROPHYLL  V  (fiq /I  )  CHLOROPHYLL 'o'  (ftg/l)  CHLOROPHYLL  'a'  (  fig  f\  ) 
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FIGURE  7 

CHLOROPHYLL  'a'  CONCENTRATION  PROJECTIONS 
LOW  FLOW  HYDROLOGY,  EFFECT  OF  DISPERSION 


FIGURE  8 

CHLOROPHYLL  V  CONCENTRATION  PROJECTIONS 
HIGH  FLOW  HYDROLOGY,  EFFECT  OF  DISPERSION 
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less  available  for  algal  uptake  in  the  lower  reservoir. 

The  effect  can  be  seen  in  the  population  decreases  in  Segments 
1,  4  and  5  starting  at  about  day  150.  As  the  upstream  algal 
concentrations  decline,  downstream  levels  are  projected 
to  increase  due  to  the  renewed  availability  of  nutrients. 

It  is  noted  that  the  algal  chlorophyll  'a*  patterns 
in  the  lower  lake  which  result  under  the  higher  dispersion 
assumption  are  quite  similar  to  those  projected  for  the 
completely-mixed  scenario.  The  situation  is  markedly  different 
for  the  lower  dispersion  assumption. 

Projections  for  the  high  flow  hydrology  are  presented 
in  Figure  8.  Chlorophyll  'a'  concentrations  are  expected 
to  range  from  the  boundary  condition  of  5  pg/1  to  peaks 
of  120  ug/1.  Overall  impoundment  concentrations  vary  depending 
upon  the  level  of  the  dispersion  coefficient  used  in  the 
computations.  The  behavior  in  the  upper  reservoir,  Segments 
2  and  3,  is  similar  to  that  for  the  low  flow  hydrology. 
Chlorophyll  'a'  concentrations  remain  low,  generally  less 
than  10  ug/1  in  that  area  until  such  time  as  flows  decrease 
to  below  500  cfs.  The  algal  population  increases  sharply 
at  that  time  due  to  increased  light  penetration  and  detention 
times.  Concentration  peaks  of  80  to  120  yg/1  are  projected 
in  that  area.  Levels  decline  thereafter  due  to  increasing 
flow  and  turbidity. 

As  under  the  low  flow  hydrology,  most  growth  occurs 
in  the  shallow  areas  of  the  impoundment,  Segments  4,  5  and 
6,  for  the  higher  dispersion  assumption.  Nutrients  brought 
into  the  lake  in  the  main  stem  segments  disperse  into  the 
bay  segments;  the  algae  produced  there  then  disperse  back 
into  the  main  stem  segments.  Peak  concentrations  in  the 
lower  lake  are  121,  117,  and  114  yg/1  for  Segments  4,  5 
and  1,  respectively.  Population  levels  peak  and  then  decline 
to  levels  that  can  be  sustained  by  nutrient  inflow  and  recycle. 
Under  the  lower  dispersion  assumption,  insufficient  nutrients 
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disperse  into  the  bay  segment.  The  shallows  are  projected 
to  exhaust  inorganic  nitrogen  concentrations  early  in  the 
year  and  not  provide  the  growth  and  biomass  exchange  seen 
in  the  higher  mixing  case.  The  result  is  a  series  of  quite 
different  and  usually  lower  annual  population  patterns  in 
the  main  body  of  the  impoundment.  The  trends  seen  in  the 
lower  reservoir  under  the  high  flow  regime  are  also  very 
similar  to  those  for  the  completely  mixed  configuration. 

2.  Dissolved  Oxygen.  The  time  variable  water  quality 
models  and  associated  segmentation  previously  discussed 
were  employed  to  generate  calculated  dissolved  oxygen  pro¬ 
files  for  the  six-segment  model  only.  The  parameters  and 
coefficients  employed  in  the  analysis  are  listed  in  Table 
3.  These  coefficients  represent  an  estimate  of  probable 
appropriate  values.  The  dissolved  oxygen  computations  were 
developed  employing  an  uncalibrated  and  unverified  model. 

The  results  of  the  analysis  should  therefore  be  viewed  as 
preliminary  estimates  of  the  anticipated  dissolved  oxygen 
levels  under  projected  conditions.  These  estimates  can, 
however,  be  employed  in  evaluation  of  the  proposed  Tennessee 
Colony  Lake  project  and  analysis  of  specific  designs  and 
operating  procedures. 

The  basic  phytoplankton  model  framework  and  values 
for  coefficients  and  transport  terms  were  employed  in  the 
dissolved  oxygen  analysis.  The  two  flow  regimes  utilized 
have  been  described  previously.  Existing  data  on  instream 
BOD  and  the  nitrogen  series  were  analyzed  resulting  in  the 
rating  curves  shown  in  Figures  4  and  5. 

The  BODg  data  and  the  line  on  the  figure  represent 
the  use  of  dissolved  oxygen  by  living  phytoplankton  during 
the  bottle  test  in  addition  to  oxidation  of  carbon  asso¬ 
ciated  with  point  source  and  other  discharges.  These  two 
sources  of  oxygen  utilization  cannot  be  separated  in  the 


TABLE  3 


COEFFICIENTS  EMPLOYED  IN  THE  LAKE  DISSOLVED  OXYGEN  ANDYS  IS 


COEFFICIENT 

2 

Dispersion  Rate  (ir.i  /day) 

Dissolved  Oxygen  Reaeration  Rate 
at  20-C  (1/day) 

BOD,  Oxidation  Rate  at  20:C 
(1/day) 

Ratio  of  Ultimate  BOD  to  BOD. 

5 

Oxygen  to  Chlorophyll  'a'  Ratio 
(mg  O^/mg  Chi) 

Phytoplankton  Respiration  Rate 
at  20°C  (1/day) 

Phytoplankton  Settling  Rate 
(m/day ) 


SYMBOL 

E 

K 

a 

K 

r 


VALUE 


c.i  or  : .  o ' : 

2.0  ?  Depth  (ft) 


C.2 


1.42 

0.133 


0.1 

0.1 


NOTE:  Coefficients  used  in  the  phytoplankton  growth  and  death  analysis 
were  developed  employing  the  calibrated  Lake  Livingston  model.  System 
geometry  and  transport  parameters  were  identical  to  those  employed  lr, 
the  phytoplankton  analysis. 


I 

I 
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standard  BOD  test  procedure.  As  the  model  converts  dying 
algae  to  an  oxygen  demand  separately  from  the  oxidation 
of  non-living  BOD^,  inputing  these  measured  BOD^  values 
would  cause  the  model  to  oxidize  the  algae  twice.  An  esti¬ 
mate  of  the  order  of  phytoplankton  contribution  to  the 
measured  river  BOD  can  be  obtained  by  consideration  of  the 
chlorophyll  levels  or  numbers  of  plankton  present  in  the 
sample.  Extensive  information  on  these  variables  is  not 
available.  However,  the  limited  available  data  suggest 
ranges  in  equivalent,  chlorophyll  between  20  and  120  ug/1 
with  average  levels  on  the  order  of  30  to  60  ig/1.  The 
calculated  equivalent  oxygen  demand  of  this  chlorophyll 
ranges  between  2  and  6  mg/1  of  five  day  BOD.  The  chloro¬ 
phyll  contribution  to  the  measured  BOD,.  data  should  be 
highest  in  the  warmer  periods  and  under  low  flow  conditions. 
In  summary,  between  one  fourth  and  one  half  of  the  measured 
BOD,,  at  Rosser  is  probably  associated  with  oxidation  of 
living  phytoplankton  biomass  in  the  BOD  bottle.  In  general, 
oxidation  of  the  plankton  biomass  will  be  upper  bounded 
by  the  BOD^  bottle  data.  If  BOD^  values  are  not  manually 
reduced  before  being  used  as  boundary  conditions,  the  dis¬ 
solved  oxygen  calculations  can  be  considered  conservative 
and  tending  to  overestimate  upstream  impacts.  An  upper 
bound  would  be  to  consider  the  bottle  BOD,,  to  be  entirely 
due  to  algae;  normal  phytoplankton  respiration  and  growth 
could  be  considered  a  lower  bound,  requiring  a  forty  percent 
reduction  in  the  measured  upstream  BOD  to  remove  algae 
effects. 

There  are  two  levels  of  waste  load  reductions  being 
considered  for  the  point  sources  in  the  Dallas-Fort  Worth 
area.  Under  1979  conditions,  it  is  anticipated  that  waste 
treatment  plants  will  be  discharging  wastewater  with 
10  mg/1  of  BOD^  and  17  mg/1  of  ammonia  nitrogen.  Further 
improvements  in  waste  treatment  under  consideration  for 
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some  future  period  would  reduce  effluent  concentrations 
to  5  mg/1  of  BOD^  and  3  mg/1  of  ammonia  nitrogen.  These 
two  projected  effluent  sets  can  be  compared  to  existing 
regionwide  effluent  quality  which  averages  approximately 
25  mg/1  of  BOD^  and  20  mg/1  of  ammonia  nitrogen.  It  is 
this  existing  level  of  treatment  which  is  reflected  in  the 
reported  BOD^.  and  ammonia  nitrogen  data  previously  illus¬ 
trated  in  Figures  4  and  5.  As  the  Trinity  River  is  domi¬ 
nated  by  sewage  treatment  plant  effluents  during  low  flow 
periods,  one  method  of  adjusting  the  two  rating  curves  to 
approximate  the  two  future  BODj.  and  NH^-N  loading  conditions 
entering  the  proposed  Tennessee  Colony  Lake  is  to  multiply 
the  curves  in  Figures  4  and  5  by  the  ratio  of  effluent  con¬ 
centrations  under  future  and  present  conditions.  These 
ratios  are  summarized  in  Table  4. 

This  procedure  is  more  appropriate  under  low  flow  condi¬ 
tions,  which  are  critical  from  the  standpoint  of  dissolved 
oxygen,  as  wastewater  contributes  a  significant  portion 
to  river  flow.  As  the  Trinity  River  flow  increases,  non¬ 
point  sources,  resuspension,  and  other  phenomena  contribute 
additional  BOD,,  and  ammonia  nitrogen  to  the  river.  Thus 
the  reductions  in  river  BOD^  and  ammonia  nitrogen  tend  to 
be  overestimated. 

The  procedure  outlined  above  includes  the  BOD^  contributed 
by  phytoplankton  oxidation  in  the  BOD  bottle  under  existing 
conditions.  Therefore,  the  calculated  dissolved  oxygen 
impacts  will  be  conservative  estimates  tending  to  yield 
low  projected  dissolved  oxygen  levels. 

Projected  dissolved  oxygen  levels  over  an  annual  period 
are  presented  in  Figures  9,  10  and  11  for  the  low  flow  regime. 
Figures  12,  13  and  14  present  comparable  information  for  the 
high  flow  regimes.  One  set  of  oxygen  profiles  was  produced 
in  response  to  upstream  BOD5  loadings  alone;  a  second  set 
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TABLE  4 

RATIOS  USED  TO  ADJUST  BODs  RATING  CURVE  AT  ROSSER 
FOR  FUTURE  LOADING  CONDITIONS 


PERIOD 

BOD,  RATIO* 

D 

Present 

1.0 

1979 

0 

O 

Long  Term 

0.20 

♦Ratio  used  to  project  curve 
to  future  conditions. 
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FIGURE  9 

DISSOLVED  OXYGEN  CONCENTRATION  PROJECTIONS 
LOW  FLOW  HYDROLOGY,  PRESENT  LOADINGS 


FIGURE  10 

DISSOLVED  OXYGEN  CONCENTRATION  PROJECTIONS 
LOW  FLOW  HYDROLOGY,  1979  LOADINGS 
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FIGURE  II 

DISSOLVED  OXYGEN  CONCENTRATION  PROJECTIONS 
LOW  FLOW  HYDROLOGY,  LONG-TERM  LOADINGS 
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FIGURE  12 

DISSOLVED  OXYGEN  CONCENTRATION  PROJECTIONS 
HIGH  FLOW  HYDROLOGY,  PRESENT  LOADINGS 
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FIGURE  13 

DISSOLVED  OXYGEN  CONCENTRATION  PROJECTIONS 
HIGH  FLOW  HYDROLOGY,  1979  LOADINGS 
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FIGURE  14 

DISSOLVED  OXYGEN  CONCENTRATION  PROJECTIONS 
HIGH  FLOW  HYDROLOGY,  LONG-TERM  LOADINGS 


of  profiles  indicates  the  combined  impact  of  the  upstream 
BOD,,  loadings  plus  the  oxygen  produced  and  used  by  phyto¬ 
plankton  growth  and  death.  The  calculations  are  performed 
assuming  vertical  uniformity  in  temperature  and  in  BOD,., 
phytoplankton  and  dissolved  oxygen  concentrations.  Strati¬ 
fication  of  the  channel  region  of  the  lake  is  possible  under 
low  flows,  however.  If  strong  vertical  density  gradients 
develop,  bottom  dissolved  oxygen  levels  may  be  gradually 
reduced,  in  time,  so  that  low  dissolved  oxygen  concentra¬ 
tions  could  be  present  in  early  fall  or  late  summer  periods. 
This  phenomenon  has  been  observed  in  Lake  Livingston. 

The  dissolved  oxygen  model  as  currently  constructed 
does  not  include  a  component  of  oxygen  demand  from  bottom 
deposits  or  from  surface  runoff  other  than  the  Trinity  River. 
A  reduction  in  the  calculated  dissolved  oxygen  level  of 
approximately  1.0  to  1.5  mg/1  could  be  associated  with 
benthal  deposits  and  other  sources  of  oxygen  demand. 

Examination  of  Figures  9  through  14  indicates  that, 
in  general,  dissolved  oxygen  levels  will  be  high  in  the 
major  portion  of  the  proposed  impoundment  represented  by 
Segments  5  and  6.  These  sections  of  the  impoundment  are 
shallow  and  vertical  stratification  should  not  cause  dis¬ 
solved  oxygen  problems.  The  sections  of  the  proposed 
Tennessee  Colony  Lake  which  are  channelized  (Segments  1, 

2  and  3)  are  projected  to  have  periods  when  average  dis¬ 
solved  oxygen  levels  drop  to  3  mg/1  under  present  loadings, 
assuming  complete  oxidation  of  incoming  phytoplankton. 
Considering  1979  loading  conditions  or  the  situation  with 
no  significant  oxidation  of  incoming  plankton,  average 
dissolved  oxygen  levels  would  reach  lows  of  5  mg/1.  Sus¬ 
tained  average  high  dissolved  oxygen  levels  above  12  mg/1 
can  be  anticipated  during  the  mid  to  late  summer.  There 
could  be  significant  diurnal  variations  in  dissolved 
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oxygen  during  these  periods.  The  calculated  dissolved 
oxygen  profiles  for  the  channel  sections  of  the  proposed 
impoundment  indicate,  for  existing  loads  as  shown  in  Figures 
9  through  14,  that  spatial  gradients  will  exist.  Lower 
dissolved  oxygen  levels  will  be  associated  with  the  upstream 
reaches.  Higher  concentrations  are  calculated  near  the 
dam.  This  situation  will  tend  to  be  offset,  to  some  extent, 
by  the  increased  depth  associated  with  the  channel  area 
near  the  dam  and  the  associated  potential  for  vertical 
stratification.  Vertical  stratification  in  the  channel 
areas  will  tend  to  result  in  low  bottom  dissolved  oxygen 
levels . 

The  information  presented  in  Figures  9  through  14  was 
developed  without  consideration  of  the  oxygen  demand  asso¬ 
ciated  with  nitrification  of  ammonia  nitrogen  in  the  pro¬ 
posed  impoundment.  Under  existing  waste  input  levels, 
substantial  concentrations  of  NH^-N  reach  the  site  of  the 
proposed  lake.  The  possibility  for  nitrification  of  this 
ammonia  exists,  and  if  it  occurs  at  significant  rates  a 
substantial  additional  reduction  in  dissolved  oxygen  can 
be  anticipated. 

There  is  no  firm  basis  for  predicting  if  nitrification 
will  occur  or  the  rate  at  which  it  will  proceed.  Nitri¬ 
fication  is  observed  in  lakes  and  streams  under  some  con¬ 
ditions.  If  nitrification  does  occur  at  high  rates,  it 
is  estimated  the  dissolved  oxygen  levels  will  be  between 
one  to  two  mg/1  in  the  channelized  portions  of  the  lake. 
These  lower  dissolved  oxygen  levels  will  extend  to  the 
region  of  the  dam.  Nitrification  generally  ceases  below 
DO  concentrations  of  1.0  to  1.5  mg/1.  The  lowering  of 
dissolved  oxygen  concentrations  due  to  nitrification  could 
be  associated  with  present  and  1979  levels  of  waste 
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discharge.  The  long-term  treatment  goals  indicated  in  Table 
4  include  nitrification  at  the  treatment  plants.  Under 
this  loading  pattern,  nitrification  would  probably  not  be 
a  significant  factor  in  the  proposed  Tennessee  Colony  Lake. 
The  impact  of  nitrification  on  the  non-channelized  portions 
of  the  impoundment  are  anticipated  to  be  minimal. 

In  the  context  of  the  assumptions  and  imperfectly  under¬ 
stood  phenomena  discussed  previously,  the  dissolved  oxygen 
situation  can  be  summarized  in  the  following  manner.  Under 
1979  waste  input  levels,  assuming  no  significant  nitrifi¬ 
cation  within  the  proposed  impoundment,  calculated  dissolved 
oxygen  levels  are  generally  in  excess  of  5  mg/1.  Periodic 
or  local  values  of  dissolved  oxygen  below  this  level  can 
be  anticipated  as  a  result  of  bottom  oxygen  demand,  vertical 
stratification  and  following  short-term  algal  blooms.  If 
nitrification  proceeds  at  high  rates  in  the  impoundment 
and  not  in  the  upstream  sections  of  the  Trinity  River, 
summer  dissolved  oxygen  levels  in  the  channelized  portions 
of  the  lake  could  drop  to  one  to  two  mg/1. 
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CHAPTER  IV 

trinity  river  water  quality 
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CHAPTER  IV 

TRINITY  RIVER  WATER  QUALITY 


A.  Purpose . 

A  second  major  concern  of  the  Phase  II  study  is  the 
impact  of  the  proposed  Tennessee  Colony  Lake  on  the  water 
quality  of  the  Trinity  River.  To  address  this  concern, 
mathematical  models  were  constructed  to  evaluate  the  water 
quality  of  the  headwaters  of  Tennessee  Colony  Lake  and  to 
evaluate  the  water  quality  impacts  in  the  river  downstream 
of  the  dam  as  a  result  of  releases  from,  the  proposed  lake. 

The  evaluation  of  the  dissolved  oxygen  profiles  in  the  river 
reaches  which  will  be  affected  by  the  backwater  from  Tennessee 
Colony  Lake  is  a  major  concern  in  this  chapter.  Another 
modeling  study  included  in  this  chapter  is  the  evaluation 
of  the  dissolved  oxygen  profiles  below  the  Tennessee  Colony 
Dam  after  releases  of  waters  of  various  quality  levels. 

B.  Milepoint  Designations. 

Several  sets  of  river  miles  of  the  Trinity  River  have 

been  cited  in  previous  studies  by  the  U.S.  Corps  of  Engineers 

and  the  Trinity  River  Authority  of  Texas  (TRA) .  Table  5 

shows  a  comparison  of  river  miles  at  selected  locations 

according  to  the  TRA  system  and  the  U.S.  Corps  of  Engineers' 

system.  Since  the  TRA  has  conducted  previous  water  quality 

f  14  1 

modeling  studies  of  the  Trinity  River,  ;  the  TRA  milepoint 
designations  are  referenced  in  this  chapter.  It  is  dif¬ 
ficult  to  maintain  a  consistent  set  of  milepoint  design¬ 
ations  since  the  length  of  the  river  has  changed  over  the 
years  as  the  result  of  oxbows  being  cut  off  and  man-made 
changes  in  channel  configuration.  Furthermore,  milepoints 
may  be  different  in  the  future  as  a  result  of  proposed 
channel  modifications. 
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C.  Previous  Modeling  Studies. 

Mathematical  models  capable  of  representing  dissolved 

oxygen  conditions  in  the  Trinity  River  upstream  of  Lake 

Livingston  were  developed  in  1974  for  the  Trinity  River 
(14) 

Authority'  The  project  report  evaluated  the  dissolved 

oxygen  conditions  in  the  Trinity  River  for  the  300  miles 
from  Fort  Worth  (Riverside)  to  below  Crockett  and  presented 
a  waste  load  allocation  analysis.  The  study  emphasizes 
the  evaluation  of  dissolved  oxygen  profiles  between  Fort 
Worth  and  the  East  Fork.  The  report  also  recommended  a 
continuing  program  to  update  water  quality  data  and  refine 
the  model  calculations. 

D.  Model  Calibration  for  1974  Conditions. 

During  the  period  of  July  11-22,  1974,  flow  in  the 
Trinity  River  reached  a  base  flow  of  45  cfs  above  the 
Riverside  reach.  During  this  time,  the  Trinity  River 
Authority  conducted  intensive  sampling  surveys  to  determine 
the  water  quality  characteristics  of  the  river,  which  were 
not  included  in  the  initial  calibration.  The  results  of 
these  sampling  efforts  provided  a  comprehensive  picture 
of  the  characteristics  of  the  Trinity  River  during  low  flow 
conditions  and  are  summarized  in  the  previous  TRA  report^14^ 
These  data  were  used  during  this  Phase  II  effort  to  extend 
the  calibration  of  the  Trinity  River  model  into  the  portion 
of  the  river  which  would  be  affected  by  the  proposed 
Tennessee  Colony  Lake.  The  model  geometry,  loadings,  and 
reaction  rates  which  were  used  are  presented  in  Tables  6 
through  8.  Comparison  of  model  results  and  measured  dis¬ 
solved  oxygen  concentrations  are  shown  in  Figure  15.  The 
model  results  provide  a  good  fit  of  the  general  trend  of 
the  data  and  can  be  considered  an  adequate  representation 
of  the  water  quality  conditions  in  the  river  during  low 
flow. 
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TABLE  5 


COMPARISON  OF  TRA  AND  U.S.  CORPS  OF  ENGINEERS 

MILEPOINTS 


LOCATION 

TRA* 

U.S.  CORPS 

OF  ENGINEERS** 

Big  Fossil  Creek 

515.3 

542.7 

Village  Creek 

506.8 

533.6 

Elm  Fork 

479.7 

505.5 

East  Fork 

435.3 

4  59.8 

Trinity  River  entering 

Seg.  3  of  Interim 

Tennessee  Colony  Lake 
(Elevation  240  ft) 

362.8 

Richland  Creek 

351.3 

372.4 

Seg.  3  to  Seg.  2  of  the 
Tennessee  Colony  Lake 

334.9 

- 

Tehuacana  Creek 

328.1 

347.2 

Tennessee  Colony  Dam 

323.5 

341.7 

Catfish  Creek 

320.7 

339.6 

Upper  Keechi  Creek 

255.3 

272.8 

Lower  Keechi  Creek 

226.1 

240.5 

♦TRA  River  Miles  are  used  in  this  report  because  previous 
modeling  works  were  based  on  TRA  River  Miles. 

♦♦Reference  2 
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FIGURE  15 

CALCULATED  AND  OBSERVED  DISSOLVED  OXYGEN  LEVELS 
IN  THE  UPPER  TRINITY  RIVER  (JULY  11-22,  1974) 


E.  Modeling  Analysis  of  the  Headwaters  of  the 
Proposed  Tennessee  Colony  Lake. 

This  section  discusses  the  application  of  a  portion 
of  the  recalibrated  river  model  to  the  headwaters  of  the 
proposed  Tennessee  Colony  Lake.  The  analysis  performed 
dealt  with  the  computed  changes  in  low  flow  dissolved  oxygen 
profiles  resulting  from  the  Trinity  River  entering  the  wider 
and  deeper  transition  section  of  the  lake.  The  impact  of 
tributary  flows  and  municipal  wastewater  treatment  plant 
effluents  entering  in  the  Dallas/Fort  Worth  area,  which 
create  loads  of  BOD,,  and  ammonia,  was  modeled.  Four  basic 
configurations  were  analyzed:  present  river  geometry  with 
and  without  instream  nitrification  of  ammonia,  and  lake- 
modified  geometry  also  with  and  without  instream  nitrifi¬ 
cation  of  ammonia.  These  four  configurations  were  analyzed 
under  several  projected  loading  patterns. 

A  detailed  discussion  of  the  mathematical  form  of  the 

model  will  not  be  presented  in  this  report.  The  model  which 

was  previously  developed  for  the  TRA  is  adequately  described 

(14) 

in  the  previous  report'  .  Instead,  only  the  particulars 
of  this  application  will  be  discussed  below,  including  model 
segmentation  and  geometry  with  the  interim  pool  of  Tennessee 
Colony  Lake,  model  kinetics,  and  the  results  of  the  dis¬ 
solved  oxygen  analyses  for  projected  future  loadings. 

1.  Segmentation.  The  most  useful  form  of  the  equation 
for  the  spatial  distribution  of  dissolved  oxygen  deficit 
and  other  water  quality  variables  in  a  river  is  derived 
assuming  constant  coefficients  for  a  specific  portion  of 
a  river.  Thus,  an  initial  step  in  river  modeling  involves 
division  of  the  river  into  reaches  of  relatively  constant 
characteristics.  The  water  quality  equation  may  then  be 
applied  separately  to  each  segment  of  the  system  to  calculate 
a  spatial  concentration  profile  within  the  respective  segment. 
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At  the  junctio  between  segments,  mass  balances  are  per¬ 
formed  to  determine  the  initial  condition  of  the  succeeding 
segment  in  the  system. 

In  summary,  a  new  segment  must  be  established  for  the 
following  reasons: 

•  at  the  location  of  any  significant  flow  discon¬ 
tinuity  in  the  system;  such  additional  flow  may 
be  in  the  form  of  a  tributary  or  wastewater 
discharge 

•  on  the  basis  of  significant  changes  in  physical 
characteristics  such  as  cross-sectional  area, 
depth,  width,  velocity,  reaeration  or  temperature 

•  when  any  of  the  reaction  rates  or  oxygen  demands 
are  known  to  vary  significantly  with  distance 

The  proposed  Tennessee  Colony  Lake  headwaters  start 
at  river  mile  362.8  when  the  conservation  pool  is  maintained 
at  an  elevation  of  240  feet.  From  this  location  down  to 
river  mile  334.9,  the  river  characteristics  led  to  a  divi¬ 
sion  of  the  headwater  reach  (Segment  3  of  the  lake  model 
presented  in  Chapter  III)  into  eight  segments.  Thus,  a 
total  of  26  segments  is  used  in  the  modified  Trinity  River 
model  for  the  evaluation  of  water  quality  in  the  headwaters 
of  Tennessee  Colony  Lake. 

2.  Geometry.  The  velocity,  depth,  width,  and  cross- 
sectional  area  for  a  specific  segment  of  a  river  can  often 
be  expressed  as  power  functions  of  flow.  The  following 
equations  are  used  to  relate  these  parameters  for  each  model 
segment : 
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Velocity: 

U 

n  NU 

-  cuQ 

(2a) 

Depth : 

H 

nh 

“  cHQH 

(2b) 

Width : 

W 

Nw 

-  cwQ 

(2c) 

Cross-sectional  Area: 

A 

na 

-  CA° 

(2d) 

in  which  Cy,  Cy,  C^,  and  C A  are  constants  of  proportionality 
and  Ny,  nh  ,  Nw,  and  NA  are  exponential  constants.  By 
reason  of  continuity,  these  coefficients  are  related  by 
the  following  expressions: 


CA 

« 

1/Cy 

(3a) 

CW 

s 

CA/CH 

(3b) 

na 

a 

1  -  Ny 

(3c) 

NW 

= 

na  -  nh 

(3d) 

Thus,  a  complete 

set 

Of 

proportionality 

constants  can 

be  determined  if  CA  or 

CU 

is 

known  and  Cw  or 

CH  is  known. 

Likewise,  the  exponential  constants  may  be  determined  if 
Na  or  Ny  is  known  and  Nw  or  Ny  is  known. 

Information  on  time  of  travel  and  river  geometry  was 
obtained  from  the  following  sources: 

•  Above  river  mile  362.8,  the  river  geometry 
developed  by  Hydroscience  in  the  previous 
modeling  effort  was  adopted 
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•  Below  river  mile  362.8,  the  average  depth  and 
width  of  the  river  with  the  lake  at  an  elevation 
of  240  feet  have  been  estimated  by  using  available 
U.S.  Army  Corps  of  Engineers'  cross-sectional  data 

The  present  exponential  relationships  for  geometry  of  the 
various  model  segments  are  shown  in  Table  9.  Revision  of 
these  relationships  can  be  made  as  additional  data  become 
available . 

3.  Kinetics.  The  reaeration  rates  (K_),  carbonaceous 

Q 

removal  and  deoxygenation  rates  (Kr  and  K^) ,  nitrogenous 
decay  rates  (Kn)  and  benthal  uptake  rates  (S)  are  parameters 
related  to  sources  and  sinks  of  the  dissolved  oxygen  in 
a  river.  These  reactions  are  described  below.  All  reaction 
rates  are  expressed  to  the  base  e. 

a.  Reaeration  Rates.  The  reaeration  coefficients 
were  calculated  using  the  O' Connor-Dobbins  formulation: 


K 


a 


(DlU) 


1/2 


H 


3/2 


(4a) 


where  is  the  oxygen  diffusivity  at  20°C  and  is  equal 
to  0.000081  ft2/hr;  U  is  average  stream  velocity  for  the 
segment;  and  H  is  the  average  depth  for  the  segment  and 
is  defined  as  the  cross-sectional  area  divided  by  the  chan¬ 
nel  width  measured  at  the  water  surface.  When  velocity 
is  expressed  in  feet  per  second  and  depth  in  feet,  this 
becomes : 


12.96  U1/2 
a  "  h3/2 


(4b) 
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to  the  base 


The  unit  of  K  in  both  formulations  is  day 

a 

e  at  20 °C. 

Where  the  Trinity  River  enters  the  headwaters  of  the 
Tennessee  Colony  Lake,  the  cross-sectional  area  and  depth 
increase  significantly.  This  would  result  in  very  low  veloc¬ 
ity  in  these  segments.  The  reaeration  rates  (K_)  for  these 

a  (15) 

segments  are  calculated  by  using  the  empirical  formula  : 


where  H  is  the  average  segment  depth  in  feet.  The  values 
of  K  in  the  headwater  reaches  of  the  lake  range  from  0.12 

3 

to  0.17/day. 

b.  Carbonaceous  Demands.  For  river  water  quality 
surveys,  carbonaceous  oxygen  demand  is  often  expressed  in 
terms  of  BOD*..  The  river  BOD  removal  rate  (Kr)  can  be  esti¬ 
mated  by  using  the  BOD,.  data  collected  at  the  various  loca¬ 
tions  of  the  river.  Based  upon  the  previous  modeling 
studies,  a  deoxygenation  rate  (Kr)  of  0.2/day  was  chosen 

for  the  model.  Because  significant  settling  occurs  at  loca¬ 
tions  below  treatment  plant  discharges,  a  slightly  higher 
total  BOD5  removal  rate  (Kf)  of  0.3/day  has  been  assumed 
in  the  model  in  the  regions  downstream  of  treatment  plants. 

c.  Nitrogenous  Demands.  The  action  of  the  nitri¬ 
fying  bacteria  can  also  consume  dissolved  oxygen  in  a  river 
through  the  oxidation  of  nitrogenous  compounds.  According 
to  the  ammonia  nitrogen  profiles  for  the  Trinity  River, 
the  decay  rate  (Kn)  of  0.0  to  0.2/day  was  selected  for  the 
modeling  of  that  parameter  for  various  portions  of  the 
river.  Generally,  dissolved  oxygen  concentrations  must 

be  above  1.5  mg/1  for  nitrification  to  occur.  This  assump¬ 
tion  affected  the  choice  of  nitrification  rates  for  this 
portion  of  the  river. 
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d.  Benthic  Demand.  Based  upon  the  previous  water 
quality  studies  of  the  Trinity,  benthic  demands  were  felt 
to  be  significant  in  certain  segments,  especially  those 
immediately  downstream  from  sewage  treatment  plants.  There¬ 
fore,  deoxygenation  due  to  benthic  demand  was  incorporated 
into  the  model.  The  values  used  in  the  model  are  based 
upon  generally  accepted  values  of  oxygen  uptake  rates  for 

municipal  sewage  sludge  in  a  river.  Values  in  the  range 

2 

of  1.0  to  4.0  gm/m  -day  are  considered  representative. 

2 

With  these  factors  in  mind,  a  value  of  1.0  to  1.5  gm/m  -day 
was  chosen  for  those  sections  with  the  major  wastewater 

2 

treatment  plant  (WWTP)  discharges  and  a  value  of  1.0  gm/m  -day 
was  selected  for  the  sections  further  downstream  of  these 
major  discharges.  These  values  include  adjustments  for 
projected  improvements  in  waste  treatment  plants  discharging 
to  the  Trinity. 

e.  Temperature  Correction  Factors.  In  the  pre¬ 
vious  discussion,  the  decay  rates,  reaeration  rates  and 
benthal  uptake  rates  were  given  to  the  base  e  at  20CC. 

In  the  various  model  runs,  appropriate  temperature  correc¬ 
tion  factors  are  applied  according  to  the  following  formulas: 


Carbonaceous  Decay: 

(Kd)T  = 

(Kd} 20 

(1.047) T_  2  0 

(5a) 

Nitrogenous  Decay: 

<Vt  = 

(V20 

(1.047)T_2° 

(5b) 

Reaeration: 

<Vt  = 

(V20 

(1.024) T_  2  0 

(5c) 

Benthal  Uptake: 

(S)T  + 

(S)  20 

(1.08)T~20 

( 5d) 

A  summer  water  temperature  of  30°C  was  used  in  the  simulations. 

4.  Point  Source  Loads.  Projections  for  point  source 
discharge  volumes  entering  the  Trinity  upstream  of  the  lake 


were  developed  for  the  years  1980  and  1990.  These  projected 
wastewater  treatment  plant  (WWTP)  discharges  and  tributary 
flows  are  shown  in  Table  10. 

By  the  year  1980,  it  is  estimated  that  10  mg/1  BOD^ , 

15  mg/1  TSS  and  10  mg/1  NH^-N  treatment  levels  will  be 
achieved  at  the  major  dischargers  to  the  Trinity  River. 

This  level  of  treatment  was  modeled  for  projected  wastewater 
volumes  under  1980  and  1990  conditions.  In  order  to  in¬ 
vestigate  some  possible  effects  of  even  higher  treatment 
levels,  the  treatment  combination  of  5  mg/1  BOD,.,  5  mg/1 
TSS  AND  3  mg/1  NH^-N  was  also  modeled  for  projected  waste- 
water  volumes  under  1990  conditions.  Table  11  summarizes 
the  effluent  concentration  and  instream  reaction  rates  for 
these  three  conditions.  Several  other  future  treatment 
and  reaction  rate  combinations  are  possible;  the  conditions 
reflected  in  Table  11  were  chosen  simply  for  the  sake  of 
example . 

5.  Dissolved  Oxygen  Analysis.  Multiple  computer  runs 
were  conducted  to  obtain  results  for  the  dissolved  oxygen 
profiles  under  the  various  loading  conditions.  Table  12 
provides,  as  an  example,  the  model  geometry  under  1980 
discharge  volumes  and  low  flow  conditions. 

Nitrification  in  a  river  is  an  uncertain  reaction  which 
may  or  may  not  occur  under  future  conditions  in  the  Trinity 
River.  As  mentioned  previously,  one  factor  which  limits 
nitrification  is  the  dissolved  oxygen  concentration. 
Furthermore,  other  phenomena  such  as  algal  uptake  or  strip¬ 
ping  may  reduce  the  ammonia  concentration  in  the  river. 

Thus,  the  nitrogen  cycle  in  a  river  such  as  the  Trinity 
is  difficult  to  evaluate.  For  this  reason,  computer  model 
results  were  obtained  for  projected  loading  conditions  with 
and  without  nitrification  of  ammonia  occurring  in  the  river. 
This  approach  hopefully  brackets  the  range  of  possibilities 
to  be  expected  in  this  portion  of  the  river. 
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TABLE  10 


PROJECTED  TRIBUTARY  FLOWS  AND  MUNICIPAL  WASTEWATER 
TREATMENT  PLANT  DISCHARGES  TO  THE 
UPPER  TRINITY  RIVER 


Upstream  (above  Fort  Worth) 

Year  1980 
(cf  s) 

50 

Year  1990 
(cf  s) 

50 

Village  Creek  WWTP 

131.5 

162.4 

TRA  Central  WWTP 

116.0 

201.0 

Elm  Fork 

22.0 

33.1 

Dallas  WWTP  (2  plants) 

278.0 

324.9 

TRA  Tenmile  Creek  WWTP 

17.0 

20.9 

East  Fork 

68.4 

125.6 

Richland  Creek 

83.0 

83.0 

UPPER  TRINITY  RIVER  MODEL  GEOMETRY  (WITH  TENNESSEE  COLONY  LAKE) 


i 
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The  DO  profiles  in  the  river  were  projected  based 
upon  a  stream  flow  which  consisted  of  upstream  flow  from 
the  headwaters  above  Dallas/Fort  Worth  and  the  waste  treat¬ 
ment  effluent  generated  in  Dallas/Fort  Worth.  An  upstream 
flow  of  50  cfs  was  used  as  the  low  flow  case.  The  pro¬ 
jected  profiles  using  a  background  of  100  cfs  and  342  cfs 
showed  a  less  severe  DO  depression  than  under  the  50  cfs 
condition.  The  results  of  the  analyses  using  50  cfs  are 
shown  in  Figures  16  through  18.  Figure  16  represents  1980 
cond itions  with  10/15/10  treatment  levels  at  all  major 
treatment  plants.  Figure  17  represents  1990  conditions 
with  10/15/10  treatment  levels  at  all  major  treatment 
plants.  The  differences  in  the  various  dissolved  oxygen 
profiles  between  these  two  years  is  minimal.  Figure  18 
represents  1990  conditions  with  5/5/3  treatment  levels 
at  all  major  plants. 

For  all  three  figures,  the  same  general  conclusion 
can  be  made:  the  backing  up  of  the  water  in  the  headwater 
region  of  Tennessee  Colony  Lake  will  apparently  cause  the 
dissolved  oxygen  to  drop  by  1.5  to  3.0  mg/1.  A  portion 
of  this  drop  would  be  caused  by  the  pollutants  remaining 
in  the  river  water  in  this  region  and  a  portion  of  this 
drop  would  be  caused  by  the  benthal  demand.  Under  present 
conditions,  in  the  area  where  the  future  headwaters  of 
Tennessee  Colony  Lake  will  be  located,  the  Trinity  River 
is  in  a  state  of  recovery  from  the  dissolved  oxygen  depres¬ 
sion  caused  by  the  pollutant  loads  in  the  Dallas/Fort  Worth 
area.  Location  of  the  lake  headwaters  in  this  region  will 
slow  down  the  velocity,  reduce  the  reaeration,  and  cause 
a  secondary  depression  of  DO,  as  shown  in  the  figures. 

For  the  assumed  conditions  without  nitrification,  the  low 
dissolved  oxygen  concentration  is  projected  to  be  between 
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FIGURE  16 

CALCULATED  DISSOLVED  OXYGEN  PROFILES 
IN  THE  UPPER  TRINITY  RIVER 
1980,  LOW  FLOW  CONDITIONS  (10/15/10) 


4  and  5  mg/1,  depending  upon  the  upstream  loading  levels. 

For  the  assumed  conditions  with  nitrification,  the  low- 
dissolved  oxygen  concentration  is  projected  to  be  between 
1  and  3  mg/1,  depending  upon  the  upstream  loading  levels. 

These  analyses  are  meant  only  to  show  the  relative 
impact  of  the  various  conditions  studied.  Numerous  other 
combinations  of  parameters  could  be  possible.  For  example, 
the  present  analysis  does  not  take  algal  activity  directly 
into  account.  Algal  activity  would  presumably  cause  wide 
diurnal  swings  in  dissolved  oxygen  in  this  portion  of  the 
river  and  several  of  the  dissolved  oxygen  profiles  shown 
could  be  raised  somewhat  on  the  daily  average  as  a  net 
result  of  algal  activity. 

F.  Modeling  Analysis  of  the  Trinity  River  Below  the 
Tennessee  Colony  Lake  Dam. 

The  second  mathematical  model  was  developed  to  evaluate 
the  water  quality  in  the  Trinity  River  brought  about  by 
the  releases  from  the  proposed  Tennessee  Colony  Lake. 

The  analysis  was  primarily  concerned  with  the  direct  effects 
of  these  releases  on  the  dissolved  oxygen  profiles  downstream 
of  the  Tennessee  Colony  Dam.  A  range  of  reasonable  ultimate 
oxygen  demand  (UOD)  and  DO  concentrations  in  the  dam  release 
water  was  used.  In  addition,  two  flow  regimes  were  studied: 
a  low  flow  of  680  cfs  and  a  moderate  flow  of  3000  cfs. 

The  predicted  effects  of  these  release  waters  of  various 
quality  may  aid  in  the  selection  of  the  release  structure 
design  and  operation. 

As  in  the  preceding  analysis,  the  model  employed  was 
developed  previously  for  the  Trinity  River  Authority^^. 

The  reader  is  referred  to  that  report  for  a  detailed  dis¬ 
cussion  of  the  model  framework.  The  following  sections 
discuss  the  model  of  the  river  downstream  of  the  dam. 

These  discussions  include  model  segmentation,  geometry, 
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kinetics,  and  evaluation  of  the  calculated  dissolved  oxygen 
profiles.  Some  of  the  technical  considerations  regarding 
the  modeling  methodology  have  been  presented  in  the  pre¬ 
vious  sections  and  these  discussions  will  not  be  repeated 
in  this  section. 

1.  Segmentation  and  Geometry.  The  model  includes 

the  portion  of  the  Trinity  River  below  the  proposed  Tennessee 
Colony  Dam  from  river  mile  323.5  to  river  mile  228.  The 
information  from  previous  modeling  work^4^  led  to  a  di¬ 
vision  of  this  part  of  the  river  into  five  segments.  A 
summary  of  the  exponents  and  coefficients  for  this  model 
is  shown  in  Table  13.  Further  geometry  for  the  two  flow 
regimes  is  given  in  Table  14. 

2.  Kinetics.  The  reaeration  rates  (K_),  carbonaceous 

cl 

removal  and  deoxygenation  rates  (Kr  and  K^) ,  nitrogenous 
decay  rates  (Kn) ,  and  benthal  uptake  rates  (S)  of  the 
Trinity  River  below  the  Tennessee  Colony  Dam  are  discussed 
in  the  following  sections. 

a.  Reaeration  Rates.  The  reaeration  rates  were 
calculated  according  to  the  O' Connor-Dobbins  formula  dis¬ 
cussed  previously  (Equation  4) .  The  calculated  rates  for 
each  segment  are  given  in  Table  15.  These  values  should 

be  reviewed  when  additional  information  on  geometry  becomes 
available. 

b.  Carbonaceous  Oxygen  Demand  and  Nitrogenous 
Oxygen  Demand.  The  combined  oxygen  demand  from  carbona¬ 
ceous  and  nitrogenous  materials  was  expressed  as  ultimate 
oxygen  demand  (UOD)  in  this  model  application.  UOD  values 
of  5  mg/1  and  15  mg/1  were  chosen  to  cover  a  reasonable 
range  of  values  possible  in  the  release  water  from  Tennessee 
Colony  Lake.  A  decay  rate  of  0.2/day  (base  e,  20°C)  was 
chosen  as  being  appropriate  for  this  portion  of  the  river. 
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c.  Benthic  Demand.  Deoxygenation  due  to  benthic 

demands  was  incorporated  into  the  model.  A  low  value  of 
2 

0.5  gm/m  -day  was  chosen  for  all  segments.  Due  to  the 
trapping  effects  of  Tennessee  Colony  Lake,  relatively  small 
organic  deposits  should  accumulate  in  the  river  bottom. 

d.  Temperature  Correction  Factors.  The  decay 
rates,  reaeration  rates,  and  benthal  uptake  rates  were 
corrected  for  temperature  according  to  the  formulas  provided 
in  the  previous  section.  A  water  temperature  of  30°C  was 
assumed  for  the  simulations. 

3.  Dissolved  Oxygen  Analysis.  Dissolved  oxygen  pro¬ 
files  in  the  river  downstream  of  the  dam  were  calculated 
for  two  flow  regimes  and  for  several  combinations  of  DO 
and  UOD  concentrations  in  the  release  water.  The  initial 
dissolved  oxygen  levels  were  chosen  as  zero  mg/1,  2  mg/1 
and  4  mg/1;  the  initial  UOD  concentrations  were  chosen  as 
5  mg/1  and  15  mg/1.  These  values  were  chosen  to  indicate 
the  expected  response  of  the  river  to  the  range  of  initial 
conditions  which  could  reasonably  be  expected  as  a  result 
of  water  releases  from  the  dam  during  low  flow  conditions. 
Figure  19  presents  the  calculated  profiles  resulting  from 
a  steady  state  release  of  680  cfs;  Figure  20  presents  com¬ 
parable  profiles  for  a  release  of  3000  cfs. 

Figure  19  shows  that  for  the  two  low  DO  cases,  the 
release  DO  concentration  controls  the  curve  and  is  the 
lowest  value  projected  in  the  river.  In  the  case  of  the 
4  mg/1  DO,  the  UOD  concentration  also  becomes  important: 
the  lower  UOD  level  appears  to  have  little  effect,  while 
the  higher  UOD  causes  a  DO  sag  to  approximately  3.5  mg/1 
at  a  distance  of  11  to  15  miles  downstream  of  the  damsite. 
Figure  20  shows  that  the  river  is  still  sensitive  to  the 
release  water  DO  and  UOD  concentrations  under  moderate  flows. 
Under  the  lower  UOD  loadings,  the  river  begins  recovering 
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FIGURE  19 

AN  EVALUATION  OF  THE  RECOVERY  CHARACTERISTICS 
OF  THE  TRINITY  RIVER  BELOW 
TENNESSEE  COLONY  DAM  (680  cfs) 
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FIGURE  20 

AN  EVALUATION  OF  THE  RECOVERY  CHARACTERISTICS 


OF  THE  TRINITY  RIVER  BELOW 
TENNESSEE  COLONY  DAM  (3000  cfs) 
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from  the  depressed  oxygen  levels  introduced  by  the  release 
waters.  Under  the  higher  UOD  loadings,  the  oxygen  concen¬ 
tration  is  further  depressed  by  the  exertion  of  the  oxygen 
demand;  the  recovery  portion  of  the  profile  does  not  begin 
until  much  farther  downstream. 

A  further  indication  of  the  slow  recovery  ability  of 
the  lower  Trinity  is  the  travel  time  required  for  a  stream 
dissolved  oxygen  level  of  5  mg/1  to  be  reached.  This  level 
is  the  most  likely  choice  of  a  stream  standard  for  this 
portion  of  the  river  in  the  future.  Table  16  summarizes 
the  approximate  times  and  distances  required  to  reach  this 
DO  level  for  the  12  combinations  of  release  flows  and  1 
loadings.  Under  low  flow  conditions,  the  standard  would 
be  met  approximately  13.5  miles  downstream  of  the  dam  in 
the  best  case,  and  54.5  miles  downstream  in  the  worst  case. 
Under  the  higher  flow,  as  the  reaeration  coefficients  are 
lowered,  the  time  required  to  recover  to  a  DO  level  of  5 
mg/1  increases,  as  does  the  equivalent  distance  downstream. 
In  the  best  case,  a  DO  level  of  5  mg/1  will  be  reached 
approximately  51.5  miles  downstream  of  the  dam;  in  the  worst 
case,  over  100  miles  downstream. 

These  results  indicate  the  sensitivity  of  water  quality 
in  the  Trinity  River  to  the  quality  of  the  Tennessee  Colony 
Lake  release  waters.  The  design  and  operation  of  the  pro¬ 
posed  release  structures  can  directly  affect  the  quality 
of  the  Trinity  River  for  a  significant  distance  downstream 
of  the  dam  by  their  control  of  the  dissolved  oxygen  concen¬ 
tration  in  the  release  waters. 
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TABLE  16 


ESTIMATED 

DISSOLVED  OXYGEN 

OF  PROPOSED 

RECOVERY 

TENNESSEE 

CHARACTERISTICS  DOWNSTREAM 

COLONY  DAM 

TEST 

RELEASE 

FLOW 

(cfs) 

TEST  RELEASE  WATER 
QUALITY 

RECOVERY  POINT' 

D.O.  1 
(mg/1) 

UOD2 3 

(mg.l) 

DISTANCE 

DOWNSTREAM  OF  DAM 
(miles ) 

TIME  OF 

TRAVEL 

(da”£) 

680 

0 

15 

54.5 

4.00 

680 

2 

15 

54.5 

4 .  o: 

680 

4 

15 

54.5 

4.c: 

680 

0 

5 

26.5 

1.7: 

680 

2 

5 

21.5 

1.42 

680 

4 

5 

13.5 

C.67 

3000 

0 

15 

>100 

>  4  . : : 

3000 

2 

15 

>100 

N  •  r-  - 

'  **  .  X  » 

.'>00 

4 

15 

>100 

>4.:: 

3000 

0 

5 

78.5 

2.59 

3000 

2 

5 

69.5 

2 . 25 

3000 

4 

5 

51.5 

1.66 

1D.O.  =  Dissolved  Oxygen 

2UOD  =  Ultimate  Oxygen  Demand  =  sum  of  long-term  carbonaceous  and 
nitrogenous  oxygen  demands 

3Point  at  which  a  D.O.  concentration  of  5  mg/1  is  achieved  following 
a  sag  or  depression  bex^.'  5  mg/1  (See  Figures  19  and  20) 
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